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Abstract 

For  the  past  five  years  there  has  been  an  ongo¬ 
ing  experimental  and  analytical  program  at  the  Uni¬ 
versity  of  Tennessee  Space  Institute  (UTSI)  to  im¬ 
prove  our  understanding  of  arcjet  physics.  A  com¬ 
putational  model  that  assumed  local  thermodynamic 
equilibrium  was  first  used  to  simulate  arcjet  thrusters 
operating  on  ammonia,  hydrogen,  and  argon.  The 
UTSI  arcjet  code  was  later  extended  to  include  a  two 
temperature,  finite  rate  kinetic  model  for  hydrogen 
plasma.  Recently,  this  code  has  been  used  to  sim¬ 
ulate  a  radiation-cooled  arcjet  (MARC  thruster)  ex¬ 
periment  and  a  water-cooled  arcjet  (TTl  thruster) 
experiment  performed  at  The  Universitat  Stuttgart 
Institut  fur  Raumfahrtsysteme.  The  results  of  these 
simulations  are  presented  along  with  a  review  of  UTSI 
arcjet  computation  code  development.  A  two-beam 
multiplexed  laser  induced  fluorescence  (LIF)  tech¬ 
nique  was  developed  at  UTSI  to  provide  detailed  mea¬ 
surements  of  arcjet  flows  near  the  nozzle  exit  plane. 
Comparison  of  detailed  flowfield  measurements  with 
predictions  of  the  computation  model  were  used  to 
provide  insight  into  the  physical  models  used  in  the 
arcjet  code.  The  method  wets  first  demonstrated  us¬ 
ing  a  small,  300  W,  water-cooled  arcjet  operated  with 
argon  propellant.  The  method  was  then  applied  to 
a  1  kW  arcjet  operated  with  hydrogen  and  nitrogen 
propellant  mixtures  using  the  Balmer  alpha  line  of 
hydrogen.  Recently,  the  method  has  been  extended 
to  use  an  excited  state  line  in  nitrogen.  The  results 
of  this  most  recent  research  are  presented. 

Introduction 

An  effort  has  been  underway  at  The  University 
of  Tennessee  Space  Institute  (UTSI)  since  1989  to 
develop  a  better  understanding  of  arcjet  physics  us¬ 
ing  computational  simulations  and  optical  diagnos¬ 
tics.  We  will  review  the  past  development  of  our  arc¬ 
jet  simulation  code  and  our  multiplexed  laser  induced 
fluorescence  (LIF)  measurements  and  then  present 
the  results  of  our  most  recent  research. 


*  B.  H.  Goethert  Professor  of  Engineering 
Science  &  Mechanics,  Member,  AIAA 
GRA,  Student  Member,  AIAA 
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The  UTSI  computational  code  was  originally 
developed  to  simulate  a  laser  sustained  plasma 
thruster  [1].  A  steady-state,  Navier-Stokes  code 
based  on  the  SIMPLE  algorithm  of  Gosman  and  Pun 
[2]  as  modified  by  Rhie  [3]  to  handle  subsonic  and 
supersonic  flows  was  extended  to  include  radiation 
transfer  and  heat  addition  [4].  The  code  proved  suc¬ 
cessful  at  predicting  the  behavior  of  laser  sustained 
plasmas  [1]  and  radio  frequency  inductive  arcs  [5],  and 
was  therefore  chosen  as  a  basis  for  further  develop¬ 
ment  to  permit  arcjet  simulations.  The  original  code 
assumed  that  the  plasma  was  at  local  thermodynamic 
equilibrium  (LTE),  which  proved  reasonable  for  laser 
sustained  plasmas  that  typically  are  operated  at  pres¬ 
sures  of  one  to  several  atmospheres.  Later,  Zerkle 
developed  a  non-equilibrium  code  for  laser  sustained 
plasmas  [6]. 

Initial  compEO^isons  of  the  UTSI  arcjet  code 
with  experiments  performed  at  NASA  Lewis  Research 
Center  [7]  using  a  10  kW  hydrogen  arcjet  showed 
that  the  code  overestimated  performance,  but  pre¬ 
dicted  correct  trends  with  propellant  flow  rate  and 
power.  Subsequent  comparisons  of  the  code  with  the 
experimental  performance  of  similar  water-cooled  and 
radiation-cooled  arcjets  revealed  the  inadequacy  of 
the  LTE  assumption  [8].  The  code  has  now  been  ex¬ 
tended  to  include  a  two  temperature,  finite  rate  ki¬ 
netic  model  for  the  plasma,  and  predictions  with  the 
new  code  are  significantly  improved.  These  new  re¬ 
sults  will  be  described  later. 

The  velocity  distribution  function  of  an  absorb¬ 
ing  atom  can  be  sampled  with  a  high  degree  of  preci¬ 
sion  using  narrow  line  tunable  dye  lasers  [9] .  For  pro¬ 
pellants  of  interest  in  arcjets,  the  only  candidates  for 
this  technique  are  excited  state  atoms,  since  the  reso¬ 
nance  absorption  lines  are  all  in  the  vacuum  ultravio¬ 
let  portion  of  the  spectrum.  This  technique  was  first 
demonstrated  in  an  arcjet  plume  using  the  Balmer 
alpha  transition  in  hydrogen  [10].  We  proposed  and 
developed  a  multiplexed  LIF  method  which  is  capable 
of  measuring  two  (or  more)  components  of  velocity 
simultaneously.  This  method  was  first  demonstrated 
using  a  300  W  arcjet  operating  with  argon  propellant 
[11].  Subsequently,  we  have  used  this  technique  to 
measure  the  exhaust  plume  velocity  in  a  1  kW  arcjet 
using  ammonia  and  hydrazine  analog  mixtures  of  hy¬ 
drogen  and  nitrogen  using  both  H-alpha  and  nitrogen 
transitions.  These  new  results  will  also  be  presented 
later. 
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2.0  The  UTSI  Arclet  Simulation  Codes 

The  general  form  of  the  UTSI  arcjet  model  is  a 
Navier-Stokes  solver  based  on  the  Semi- Implicit  Pres¬ 
sure  Linked  Equation  (SIMPLE)  algorithm  of  Gos- 
man  and  Pun  [2],  which  was  modified  by  Rhie  to  han¬ 
dle  subsonic  and  supersonic  fiows  [3].  The  capabili¬ 
ties  for  heat  addition  and  radiation  heat  transfer  were 
later  added  by  Jeng  and  Keefer  [4],  The  UTSI  arcjet 
model  has  been  extended  to  remove  the  assumptions 
of  equilibrium  chemical  composition  and  LTE.  The 
latter  assumption  has  been  replaced  with  a  two  tem¬ 
perature  model  for  the  electrons  and  the  heavy  gas. 
With  a  few  significant  exceptions,  the  current  formu¬ 
lation  for  the  non-equilibrium  model  was  described 
in  a  paper  presented  at  the  23rd  International  Elec¬ 
tric  Propulsion  Conference  [12],  Since  that  work  was 
present^,  the  model  has  also  been  modified  to  re¬ 
move  the  constraints  of  specified  power  dissipation 
and  specified  current  distribution  on  the  anode. 

2.1  The  Equilibrium  Arcjet  Simulation  Code 

Numerical  modeling  of  the  fluid  dynamics  in  an 
arcjet  is  a  complicated  task.  Three  momentum  equa¬ 
tions,  the  energy  equation,  and  the  continuity  equa¬ 
tion  must  be  solved  simultaneously  with  three  elec¬ 
tromagnetic  equations.  The  UTSI  arcjet  model  teikes 
advantage  of  several  assumptions  and  simplifications 
to  make  the  problem  manageable. 

The  UTSI  code  [13]  uses  the  fluid  equations  for 
steady,  axially  symmetric,  compressible  flow.  In  the 
equilibrium  model,  local  thermodynamic  and  chem¬ 
ical  equilibrium  is  assumed  everywhere.  This  as¬ 
sumption  implies  that  heavy  particle  and  electron 
temperatures  are  equal  and  that  all  thermodynamic 
and  transport  properties,  including  species  concen¬ 
trations,  can  be  defined  by  two  state  variables.  The 
flow  is  assumed  to  be  laminar.  The  electromagnetic 
equations  are  also  steady  in  time  and  axially  symmet¬ 
ric.  The  axial  and  radial  components  of  the  magnetic 
field  aire  assumed  to  be  zero. 

One  exception  to  the  assumption  of  local  ther¬ 
modynamic  equilibrium  (LTE)  is  made  to  account 
for  non-equilibrium  effects  of  electrical  conductivity. 
Changes  to  the  electrical  conductivity  can  greatly  in¬ 
fluence  the  character  of  the  flow  because  this  property 
controls  the  current  distribution  and  the  electrical 
power  absorption  in  the  flow.  The  assumption  of  LTE 
caiises  the  conductivity  to  be  a  function  of  pressure 
and  static  enthalpy.  However,  in  certain  areas  of  the 
flow,  such  as  in  the  downstream  expansion  and  near 
the  anode  wall,  the  gas  is  not  in  equilibrium.  In  par- 
ticular,  assuming  LTE  in  the  boundary  layer  near  the 
anode  seriously  changes  the  physics  of  the  model.  Un¬ 
der  true  LTE  conditions,  the  relatively  cold  boundary 
layer  would  have  an  extremely  low  conductivity  that 
would  raise  the  voltage  at  the  wall  causing  the  charac¬ 
ter  of  the  arc  reattachment  to  be  unrealistic.  Under 
actual  conditions,  the  boundary  layer  at  the  anode 


has  enhanced  electron  concentration  due  to  chemical 
and  thermodynamic  non-equilibrium  in  the  boundary 
layer.  Compensation  for  the  non-equilibrium  effects 
is  achieved  by  artificially  enhancing  the  electrical  con¬ 
ductivity  at  the  anode.  Near  the  wall,  the  conductiv¬ 
ity  is  not  allowed  to  decrease  at  low  temperatures  as 
it  would  under  an  LTE  calculation.  Figure  1  shows 
this  effect  graphically. 


Figure  1.  Electrical  conductivity  modification  used  in  the 
LTE  version  of  the  simulation  code. 


2.1.1  Governing  Equations 

The  model  solves  for  eight  dependent  variables 
from  eight  coupled  partial  differential  equations.  The 
dependent  variables  are  axial  velocity  (n),  radial  ve¬ 
locity  {v),  azimuthal  velocity  (tt;),  static  enthalpy 
(iif),  pressure  (p),  the  azimuthal  component  of  mag¬ 
netic  field  {Bq),  and  the  axial  and  radial  current  den¬ 
sities  (Jx)  and  (Jr).  The  electromagnetic  equations 
are  coupled  to  the  fluid  equations  through  Ohmic 
heating  and  Lorentz  body  forces.  The  following  sym¬ 
bols  are  used  in  the  governing  equations: 

X  axial  coordinate 

r  radial  coordinate 

p  gas  density 

/i  viscosity 

Cp  specific  heat  at  constant  pressure 

k  thermal  conductivity 

a  electrical  conductivity 

Qi  Ohmic  heating 

Or  radiation  loss 

$  viscous  dissipation 

V  vector  velocity 

The  velocity  components  are  obtained  from  con¬ 
servation  of  momentum  in  three  directions.  The  mo- 
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mentum  equation  in  vector  notation  is: 


J  =  (t[E+ (v  XB)]  (2.10) 


V-pv-v  =  V-pVv— Vp— V  X  (pV  xv)+JxB 

(2.1) 

The  static  enthalpy  is  obtained  from  the  energy  equa¬ 
tion 


k 

V-pvH  =  V — VH+^+qi-qr+v'^p  (2.2) 

Cp 

which  includes  the  viscous  dissipation  term.  For  ax- 
isymmetric  flow  the  viscous  dissipation  is, 


2  /  dv  V  ^ 

3  V  ^  r  / 

(2.3) 

The  Ohmic  heating  term  is  given  by, 

+  .  (2.4) 

(T 

The  radiation  transport  in  the  energy  equation  is 
modeled  by  assuming  optically  thin  and  optically 
thick  components  as  us^  by  Jeng  and  Keefer  [4]. 
The  optically  thick  radiation  is  modeled  using  the 
Rosseland  [14]  approximation  which  adds  a  term  to 
the  local  thermal  conductivity  of  the  gas.  The  opti¬ 
cally  thin  radiation  is  treated  as  energy  lost  from  the 
system.  The  pressure  is  obtained  from  the  continuity 
equation. 

V  •  pv  =  0  (2.5) 

The  electromagnetic  properties  of  the  flow  are  de¬ 
scribed  by  Maxwell’s  equations  and  Ohm’s  law  using 
the  MHD  approximation  which  neglects  displacement 
current  and  assumes  that  the  fluid  is  electrically  neu¬ 
tral: 

V  X  B  =  poJ  (2.6) 

V  X  E  =  0  (2.7) 

V  •  B  =  0  (2.8) 

V  •  E  =  0  (2.9) 


where  /jLq  is  the  permeability  of  free  space,  O  is  the 
electrical  conductivity  and  B,  E,  v,  and  J  are  the 
vector  quantities  of  magnetic  field,  electric  field,  ve¬ 
locity,  and  current  density.  These  equations  are  com¬ 
bined  and  simplified  by  assuming  axial  symmetry  and 
by  assuming  that  the  radial  and  axial  components  of 
magnetic  field  are  zero.  The  resulting  equation  for 
azimuthal  magnetic  field  written  in  terms  of  the  com¬ 
bined  variable  (3  =  tBq  is. 


and  the  components  for  current  density  are  written. 


/zor  dr 
-1  dp 
flQT  dx 


(2.12) 

(2.13) 


The  governing  Eqs.  2.1  through  2.13  are  solved  in  a 
transformed  coordinate  system,  which  is  a  mapping 
of  the  true  nonorthogonal  and  nonuniform  grid  to  a 
grid  of  unit  squares.  Nonlinear  terms  are  linearized  by 
including  them  in  source  terms  [15].  Equation  2.11, 
which  describes  the  magnetic  field,  is  lineEir  in  /5.  It  is 
solved  by  direct  matrix  inversion.  A  solution  for  (3  is 
typically  calculated  every  20  iterations.  The  solution 
is  relaxed  into  the  previous  solution  with  Eq.  2.14. 
The  fluid  equations  for  IZ,  V,  tU,  if,  and  p  are  solved 
sequentially  using  updated  values  of  the  dependent 
variables  when  they  are  available  and  the  old  values 
when  updated  ones  have  not  yet  been  calculated.  As 
the  new  values  of  the  dependent  variables  are  calcu¬ 
lated,  they  also  are  relaxed  into  the  previous  solution 
according  to  Eq.  2.14, 


—  (1  ““  fr)  ^o\d  +  /r^new  •  (2.14) 

In  Eq.  2,14,  ({>  is  any  dependent  variable 
if,p,  or  0)  and  relaxation  factor  between 

zero  and  one.  A  typical  value  for  fr  is  0.2.  Solu¬ 
tion  iterations  continue  until  the  residual  errors  of 
the  equations  have  been  reduced  to  reasonably  small 
values.  Upwind  differencing  for  convective  terms  and 
fourth  order  pressure  smoothing  maintain  stability  of 
the  dependent  variables  before  a  final  solution  is  de¬ 
termined  [3]. 
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2.1.2  Boundary  Conditions 

A  no-slip  velocity  boundary  condition  is  in  ef¬ 
fect  at  all  walk.  The  wall  temperatures  are  specified  a 
priori  at  each  axial  location,  and  the  normal  pressure 
gradient  on  the  walk  k  zero.  The  radial  components 
of  U  and  V  and  the  derivatives  of  U,  H,  and  p  are 
set  to  zero  along  the  symmetry  line  at  the  center  of 
the  flow.  The  axial  derivatives  of  all  dependent  vari¬ 
ables  are  set  to  zero  at  the  nozzle  exit.  Velocity  and 
enthalpy  are  specified  at  the  inlet  based  on  the  mass 
flow,  inlet  temperature,  and  inlet  swirl  angle,  which 
are  input  parameters  to  the  program.  The  boundary 
condition  for  the  magnetic  field  k  specified  with  the 
combined  variable  /3.  The  boundary  condition  on  /3  k 
determined  from  the  axial  component  of  the  current 
by  integration  of  Eq.  2.12  with  respect  to  r  from  the 
center  of  the  flow  to  the  wall  at  each  axial  location. 
Thk  yields  the  equation 

/^wall  “  (-S0r)wall  =  /^O  /  Jx^dv  = 

(2.15) 

where  I  k  the  total  current  inside  the  boundary  of 
interest  at  a  particular  eocial  location.  Therefore,  (3 
can  be  specified  on  the  cathode  and  anode  walk  if 
the  current  /  k  known  at  each  axial  location.  The 
current  in  the  arcjet  originates  at  the  cathode  tip, 
and  k  assumed  to  flow  axially  and  uniformly  into  the 
gas  over  an  area  of  about  one  percent  of  that  of  the 
constrictor.  The  current  flow  into  the  nozzle  wall  k 
assumed  to  be  linearly  dktributed  in  the  expansion 
section  of  the  nozzle.  Current  k  not  allowed  to  ex¬ 
tend  past  the  nozzle  exit  plane.  These  assumptions 
are  based  on  observations  of  arcjets  in  operation  [16]. 
Because  a  value  for  I  can  be  specified  for  all  axial 
locations,  the  magnetic  field  can  be  calculated  on  all 
boundaries. 

Recently,  the  assumption  of  a  linear  distribution 
of  current  on  the  anode  has  been  relaxed  by  changing 
the  boundary  condition  on  /3  from  a  specified  value 
to  a  zero  normal  gradient,  starting  just  downstream 
of  the  end  of  the  constrictor.  Thk  has  the  effect  of 
putting  a  constant  voltage  boundary  condition  on  the 
nozzle  and  allows  the  current  dktribution  to  adjust 
depending  on  the  plasma  conductivity.  Thk  k  a  more 
realktic  representation  of  the  boundary  physics,  and 
results  in  a  significant  change  in  the  current  dktribu¬ 
tion  on  the  anode. 

A  limitation  k  placed  on  the  total  electrical 
power  absorbed  by  the  flow.  The  model  calculates 
local  power  absorption  from  local  current  and  con¬ 
ductivity.  The  model  assumes  LTE,  which  results  in 
values  of  electrical  conductivity  in  the  flow  which  may 
differ  considerably  from  those  in  a  non-equilibrium 
state.  The  conductivity  dictates  the  power  absorp¬ 
tion  which,  in  turn,  affects  the  entire  flow.  To  allow 
for  the  general  uncertainty  in  the  electrical  conductiv¬ 
ity,  the  sum  of  calculated  absorbed  power  from  each 


computational  location  in  the  flow  k  forced  to  equal 
a  predetermined  total  power.  Thk  k  accomplkhed 
by  multiplying  the  power  absorbed  at  each  computa¬ 
tional  location  by  the  ratio  of  the  specified  total  power 
input  to  the  absorbed  power  summed  throughout  the 
flow.  Thk  correction  insures  that  the  correct  total 
power  k  deposited  into  the  flow.  The  correction  was 
first  added  to  a  precursor  to  the  UTSI  arcjet  code,  a 
code  developed  to  simulate  inductive  radio  frequency 
arc  discharges,  to  aid  in  stability  [5|. 

The  constraint  on  total  power  has  also  been  re¬ 
laxed  for  our  most  recent  calculations.  Since  power 
supplies  usually  regulate  current,  total  power  k  an 
experimental  result,  and  it  k  more  realistic  to  allow 
the  calculated  results  to  predict  the  total  power.  For 
these  new  calculations,  the  total  current  k  specified, 
and  the  total  power  k  determined  by  the  integral  of 
Qi  over  the  volume  of  the  arcjet. 

2.1.3  Equilibrium  Thermodynamic  and  Transport 
Properties 

The  thermodynamic  and  transport  properties 
used  by  the  code  are  stored  in  lookup  tables  that  the 
program  can  access  during  operation.  Temperature 
k  stored  as  a  function  of  static  enthalpy  for  a  range 
of  pressures.  The  transport  properties  vkcosity  (/i), 
thermal  conductivity  (A:),  and  electrical  conductiv¬ 
ity  (a)  are  stored  as  functions  of  temperature  and 
pressure  [15]. 

The  ammonia  property  tables  were  calculated 
by  R.  J.  Zollweg  [17]  using  collision  cross  section 
data  for  a  hydrogen  and  nitrogen  mixture.  The  cal¬ 
culations  were  mauie  for  a  non-ideal  plasma  in  local 
thermodynamic  equilibrium.  The  electrons  in  a  non¬ 
ideal  plasma  have  kinetic  energies  of  the  order  of  the 
Coulomb  potential,  so  that  only  a  few  electrons  are 
inside  a  Debye  sphere  [18].  Zollweg  writes  that  non¬ 
ideal  plasmas  occur  “at  relatively  low  temperatures 
but  relatively  high  electron  densities.”  These  condi¬ 
tions  can  exkt  over  a  significant  portion  of  the  arcjet 
flowfield  [18]. 

Local  thermodynamic  equilibrium  (LTE)  and 
chemical  equilibrium  are  assumed  in  the  hydrogen 
equilibrium  code.  Therefore,  temperature  and  species 
concentrations  are  functions  of  static  enthalpy  and 
pressure  only.  The  temperature  k  tabulated  as  a 
function  of  enthalpy  for  a  range  of  pressures  and  k 
obtained  by  table  lookup  as  needed.  The  transport 
properties:  viscosity,  thermal  conductivity,  and  elec¬ 
trical  conductivity,  are  tabulated  as  functions  of  tem¬ 
perature  and  pressure.  The  program  used  to  generate 
these  properties  [19]  uses  tabulated  collision  cross  sec¬ 
tion  data  for  the  H2^  H'^^e  system  and  outputs 

the  internal  thermal  conductivity,  electrical  conduc¬ 
tivity,  vkcosity,  and  the  binary  diffusion  coefficients. 
Ambipolar  diffusion  is  assumed  for  the  electrons  since 
the  model  assumes  no  charge  separation. 

The  radiation  from  hydrogen  as  a  function  of 
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temperature  and  pressure  at  chemical  and  thermody¬ 
namic  equilibrium  were  calculated  using  the  methods 
described  by  Griem  [20].  Energy  transport  by  ra¬ 
diation  is  modeled  as  an  optically  thin  component 
which  results  in  a  direct  energy  loss  from  the  system, 
and  an  optically  thick  component  which  is  used  to 
define  an  equivalent  thermal  conductivity  using  the 
Rosseland  approximation  [14],  An  effective  thermal 
conductivity  is  calculated  as  the  sum  of  the  internal 
conductivity,  the  conductivity  from  the  Rosseland  ap¬ 
proximation,  and  the  energy  transport  resulting  from 
the  diffusion  of  species  [21]. 

2.2  The  Non-equilibrium  Arcjet  Simulation  Code 

The  LTE  model  for  the  plasma  is  replaced  by  a 
model  which  provides  separate  temperatures  for  the 
electrons  and  heavy  particles  and  finite-rate  chemical 
kinetics.  The  form  of  the  momentum  and  continuity 
equations  are  the  same  for  both  the  equilibrium  and 
non-equilibrium  models.  The  density  used  in  these 
equations  is  modified  to  account  for  the  two  tem¬ 
perature  plasma.  Starting  with  the  definition  of  the 
partial  pressure  for  any  species  (j)i  =  UikTi),  the 
mixture  density  (p  =  Y^niTTli),  and  the  mbcture 
pressure  (p  =  mixture  density  may  be 

written: 

p  =  pMfR[Th  +  {Te-Th)xe]  ,  (2.16) 

where  M  is  the  mixture  molecular  weight,  R  is  the 
gas  constant,  and  Xg  is  the  electron  mole  fraction. 

2.2.1  Mixture  Enthalpy 

In  equilibrium  flow,  the  energy  transported  by 
diffusion  of  species  may  be  combined  with  the  term 
containing  fcVT,  since  the  species  gradients  may  be 
written  in  terms  of  the  temperature  gradient.  For 
flow  that  is  not  in  equilibrium,  these  terms  must  be 
considered  separately.  In  this  case,  the  equation  for 
mixture  enthalpy  becomes: 

V  -  pvH 

=  V . 

+  (j>D-  k/cpj  E  /liVoii  -  kaeCpJCp^VTe] 
+  $  +  -  gr  +  v- Vp  (2.17) 

where  Cp^  is  the  sum  of  Cp^ai  over  the  heavy  species, 
and  subscript  6  represents  electrons. 

2.2.2  Electron  Energy 

The  assumption  of  LTE  made  in  our  previ¬ 
ous  work  is  relaxed  in  the  non-equilibrium  model  by 
assuming  the  electron  gas  and  the  heavy  particles 


have  separate  temperatures,  both  with  a  Maxwell- 
Boltzmann  distribution.  This  requires  a  separate  en¬ 
ergy  equation  for  the  electron  gas.  In  this  equation 
the  electrical  Ohmic  heating  is  balanced  by  electron 
convection,  conduction,  diffusion,  and  energy  trans¬ 
fer  to  the  heavy  particles  by  collisions.  The  form  of 
the  equation  and  the  included  terms  are  taken  from 
Park  [22], 

V  •  pvaehe  -  V  •  [fceVTe  +  /icDVae] 

-vVpe  =  5'  (2.18) 

where  S  =  Qce  +  Qch  +  Qi  —  Qr- 

The  term  describing  the  work  done  by  the 
electron-pressure  gradients  (v*  Vpe)  is  obtained  using 
the  assumption  that  the  electrons  and  heavy  species 
share  the  same  flow  velocity.  The  electron  partial 
pressure  is, 

Pe  =  pnxJiTh  +  (Te  -  Th)x,]  (2.19) 

where  Xe  is  the  mole  fraction  of  electrons. 

The  term  Qce  represents  the  energy  transfer  be¬ 
tween  the  heavy  gas  and  the  electron  gas  resulting 
from  collisions  of  electrons  with  the  heavy  species, 

q^  =  y^2nJ^n\k(T-T,)  .  (2.20) 

^  rrii  Z 

The  collision  frequency  (l^i)  is  evaluated  from  the 
temperature  dependent  collision  cross-sections  for 
electron-heavy  collisions  evaluated  at  the  electron 
temperature  using  the  transport  properties  program 
described  in  Section  2.2.5. 

The  term  qch  represents  the  energy  transport 
to  the  electron  gas  resulting  from  the  participation  of 
electrons  in  chemical  reactions.  It  is  assumed  that  in 
reactions  where  an  electron  is  the  third  body  the  dis¬ 
sociation  or  ionization  energy  comes  from  the  electron 
gas,  lowering  the  electron  temperature,  and  the  en¬ 
ergy  of  recombination  goes  to  the  electron  gas,  raising 
the  electron  temperature. 

The  energy  input  to  the  system  ((Ji)  and  the 
radiation  from  the  system  (^r)  are  evaluated  in  the 
same  way  as  they  are  in  the  equilibrium  model. 

2.2.3  Species 

Three  species  equations  were  added  to  the  sys¬ 
tem  of  equations  to  allow  the  calculation  of  flows 
where  dissociation  of  H2  and  ionization  of  H  are  not 
in  equilibrium.  The  mixture  is  assumed  to  consist  of 
four  species:  molecular  hydrogen  {H2),  atomic  hy¬ 
drogen  (H),  ionized  hydrogen  and  electrons 

(e) .  With  the  assumption  of  electrical  neutrality,  the 
electron  and  ion  concentrations  are  equal,  and  one 
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species  equation  can  be  eliminated.  The  equation  for 
conservation  of  a  species  may  be  written  as, 

V  ^  /Tva  =  V  •  pDVa  +  5a  ,  (2.21) 

where  So.  is  the  production  of  the  species  by  chemi¬ 
cal  reaction.  A  residence  time  is  evaluated  from  the 
mass  in  a  computational  cell  divided  by  flux  into  the 
cell,  and  the  species  production  equations  are  solved 
over  this  residence  time  step  using  a  linearized  set  of 
equations  for  the  four  species  [23].  The  source  term 
then  becomes  the  change  in  the  species  concentration 
divided  by  the  time  for  the  step. 

2.2.4  Boundary  Conditions 

Additional  boundary  conditions  are  required  for 
the  electron  temperature  and  the  species  equations 
for  the  nonequilibrium  flow  model.  The  walls  are  as¬ 
sumed  to  be  catalytic  for  electron-ion  and  H  atom  re¬ 
combination,  have  a  no-slip  boundary  condition  and  a 
specified  temperature  distribution.  As  a  consequence, 
the  gas  is  in  equilibrium  at  the  wall  temperature. 
Transport  properties  at  the  wall  are  evaluated  at  the 
conditions  of  the  first  computational  point  inside  the 
wall.  The  gas  is  assumed  to  be  in  equilibrium  at  the 
inlet,  although  a  minimum  mole  fraction  of  10“^°  is 
set  for  any  species  throughout  the  flowfield.  The  two 
temperatures  are  constrained  to  a  minimum  value  of 
100  K. 

2.2.5  Transport  Properties  for  Non-equilibrium 

Hydrogen 

Transport  properties  in  a  two  temperature 
model  with  non-equilibrium  species  concentrations 
depend  on  the  two  temperatures  and  the  number  den¬ 
sities  of  all  the  species.  For  non-equilibrium  hydrogen 
the  properties  are  a  function  of  five  variables  repre¬ 
senting  the  two  temperatures,  the  total  number  den¬ 
sity,  and  the  degree  of  dissociation  and  ionization. 
Both  a  five  dimensional  table  lookup  and  direct  cal¬ 
culation  at  each  field  point  are  impractical,  so  we  de¬ 
cided  to  combine  curve  fitting,  table  lookup,  and  some 
simplifying  assumptions  to  give  a  procedure  that  is 
both  accurate  and  reasonably  fast.  One  requirement 
is  that  the  properties  vary  smoothly  over  a  wide  range 
of  variables  to  avoid  instabilities  while  the  solution  is 
converging. 

The  program  we  are  using  to  obtain  the  trans¬ 
port  properties  [19]  uses  theoretical  and  experimental 
collision  cross-section  data  to  obtain  collision  inte¬ 
grals  from  which  the  transport  properties  are  calcu¬ 
lated.  It  was  found  that  the  transport  properties  at 
constant  pressure  and  species  mole  fraction  could  be 
fitted  to  an  equation  of  the  form  Z  =  q!T^  where  T' 
is  the  electron  temperature  for  electron  thermal  con¬ 
ductivity  {ke)  and  electrical  conductivity  (a),  and  T 
is  the  heavy  gas  temperature  for  heavy  gas  thermal 
conductivity  {kh),  and  viscosity  (fi). 


The  thermal  conductivity  of  the  electron  gas 
used  in  the  electron  energy  equation  is  taken  to  be 
Jtg.  Park  [22]  shows  that  this  value  may  be  too  high, 
since  part  of  the  energy  transfer  implicit  in  this  coef¬ 
ficient  is  transferred  to  the  heavy  particles. 

For  ke  and  (7  the  coefficients  a  and  n  are  func¬ 
tions  of  pressure  and  electron  concentration,  but  were 
found  to  be  relatively  insensitive  to  the  degree  of  dis¬ 
sociation  of  H2-  Therefore,  the  H2IH  ratio  is  not 
accounted  for  in  the  fit  coefficients.  The  properties 
kh  and  /i  are  calculated  for  each  species  and  are  com¬ 
bined  for  the  mixture.  For  the  electron-ion  pair  the 
fit  coefficients  were  found  to  be  functions  of  pressure, 
while  for  H2  and  H  they  are  constants.  The  best 
way  found  to  get  the  mixture  properties  is  to  sum  the 
product  of  the  mole  fraction  and  the  species  property. 
For  the  present  case,  this  procedure  gave  results  much 
closer  to  those  calculated  for  the  actual  mixture  than 
did  the  Wilke  mixture  rule  [24]. 

The  procedure  described  above  was  checked  by 
comparing  its  results  with  those  for  equilibrium  mix¬ 
tures  using  the  transport  property  program.  The 
agreement  is  within  a  few  percent  for  all  the  cases 
checked. 

The  binary  diffusion  coefficients  {'Dij)  are  equal 
to  cT^/p  where  C  is  a  function  of  p  for  'Dh,H2 
a  constant  for  coeffi¬ 

cients  involving  ions  are  mutiplied  by  (1  +  Te/Th) 
to  account  for  ambipolar  diffusion.  In  these  flows 
there  are  few  locations  where  H2  and  the  electron- 
ion  pair  coexist.  A  single  diffusion  coefficient  for  the 
mixture  was  used,  rather  than  the  much  more  compli¬ 
cated,  but  correct,  method  of  calculating  multicom¬ 
ponent  diffusivities.  To  provide  a  smooth  transition 
across  regions  where  all  three  components  coexist,  a 
heurestic  mixture  rule  is  used, 

D  =  (0^10:2^12  +  0^10^3^13  +  0L20tzT^2zl 
(aiQ:2  +  0:10:3  +  0:20:3)  .  (2.22) 

This  gives  the  correct  answer  for  a  mixture  of  any  two 
components  and  provides  a  smooth  transition  where 
three  component  mixtures  exist. 

3.0  Multiplexed  Laser  Induced  Fluorescence 
Velocimetry 

A  two-beam  multiplexed  laser  induced  fluores¬ 
cence  (LIF)  technique  was  developed  to  provide  de¬ 
tailed  measurements  of  arcjet  flows  near  the  nozzle 
exit  plane.  Comparison  of  detailed  flowfield  measure¬ 
ments  with  predictions  of  the  computational  model 
were  used  to  provide  insight  into  the  physical  models 
used  in  the  code  [25].  The  method  was  first  demon¬ 
strated  using  a  small,  300  W,  water-cooled  arcjet  op¬ 
erated  with  argon  propellant  [11].  The  method  was 
then  applied  to  a  1  kW  arcjet  operated  with  hydrogen 
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Balmer  alpha  line  of  hydrogen.  Recently,  the  method 
has  been  extended  to  use  an  excited  state  line  in  ni¬ 
trogen  [26]. 

The  Doppler  principle  for  measuring  flow  veloc¬ 
ities  with  LIF  has  been  known  for  some  time.  Zim¬ 
merman  and  Miles  [27]  reported  an  LIF  technique  to 
measure  the  velocity  in  a  hypersonic  wind  tunnel  in 
1980.  A  narrow  line,  tunable  dye  laser  is  used  to 
excite  an  atom  in  the  flow.  Fluorescence  from  the  ex¬ 
cited  atom  is  detected  as  the  laser  is  tuned  across  the 
absorption  line  profile.  The  wavelength  at  which  the 
line  is  excited  depends,  through  the  Doppler  effect, 
on  the  velocity  of  the  atom.  If  the  spectral  width  of 
the  laser  is  small  compared  to  the  absorption  Doppler 
profile  of  the  atom,  then  the  measured  line  profile  is 
a  measure  of  the  velocity  distribution  function  of  the 
atoms.  Thus,  in  principle,  mean  velocity  along  the 
laser  beam  direction  and  temperature  can  be  deter¬ 
mined  from  the  measured  Doppler  line  profile. 

In  the  multiplexed  LIF  technique  the  l^r  beam 
is  divided  into  two  beams  which  are  recombined  with 
an  angle  between  them.  Each  beam  is  modulated  at 
a  different  frequency  and  phase-lock  detection  is  used 
to  separate  the  fluorescence  induced  from  each  beam. 
The  velocity  along  each  of  the  beam  directions  can  be 
combined  to  provide  a  simultaneous  vector  measure 
of  velocity.  One  difficulty  with  this  technique  is  that 
an  accurate,  unshifted  reference  wavelength  must  be 
measured  to  obtain  an  accurate  absolute  velocity.  We 
solved  this  problem  by  using  an  optogalvanic  cell  to 
pj-Qyide  a  continuous,  simultaneous  reference  using 
the  same  atomic  transition  [28].  A  schematic  of  a 
typical  experimental  arrangement  is  shown  in  Fig.  2. 


Figure  2.  Schematic  of  the  complete  LIF  setup:  M  —  mir¬ 
ror,  BS  =  beam  splitter,  DFC  =  dual  frequency  chopper, 
HC  =  hollow  cathode  lamp,  LI  =  focusing  lens,  L2  =  col¬ 
lection  lens,  PD  =  photodiode  (with  pinhole  aperture), 
_  arcjet,  and  LOCK  —  phase-iocked  amplifier. 


3.1.1  Descritition  of  the  Multiplexed  LIF 

Technique 

A  detailed  theoretical  analysis  of  the  multi¬ 
plexed  LIF  technique  is  given  by  Ruyten  and  Keefer 
[11].  Here  we  will  present  only  the  elements  essential 
to  understanding  the  basic  process. 

It  is  assumed  that  the  velocity  distribution  of 
the  diagnostic  species  in  the  arcjet  plume  is  described 
by  /(Ux)  Uy)  .  where  Vx  and  Vy  are  the  velocity  com¬ 
ponents  parallel  and  perpendicular  to  the  arcjet  axis. 
The  mean  velocity  components  <  Vx  >  and  <  Uy  > 
are  given  by  the  first  moments  of  the  distribution 
function. 


oo  oo 
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The  Doppler  shift  in  velocity,  VD  is  given  by, 


VD  =  —  ^o)  (2-^) 


where  A  is  the  wavelength  of  the  transition,  Ui,  is 
the  laser  frequency  and  Uq  is  the  center  frequency  of 
the  unshifted  atomic  transition.  If  the  homogenous 
linewidth  of  the  absorption  process  (laser  linewidth 
and  homogeneous  atomic  linewidth)  is  small 
pared  to  the  Doppler  linewidth,  then  the  LIF  signal 
produced  by  the  beam  at  angle  (f>  to  the  centerline  is 
given  by, 


S4vD)  =  C4,f  [  f{vx,vy)dvxdvy  (3.3) 

J  Jband 


where  is  an  experiment-dependent  constant.  The 
“band”  of  integration  is  given  by  the  Doppler  reso¬ 
nance  condition. 


\vx  COS VySin(j)  -  vd\  ^2^^ 

where  Al/  is  the  small  homogeneous  linewidth  of  the 
excitation  process.  Since  the  homogeneous  linewidth 
is  much  smaller  than  the  Doppler  width,  the  signal  is 
essentially  proportional  to  the  component  of  the  ve¬ 
locity  distribution  function  in  the  direction  (p.  Thus, 
it  is  possible  using  Eq.  3.1  to  determine  the  mean 
velocity  <  >  in  the  directions  <f)  for  both  beams. 

In  practice,  one  beam  is  directed  perpendicular 
to  the  arcjet  axis  and  the  second  beam  is  direct^ 
along  an  angle  (f)  chosen  to  provide  a  Doppler  shift 
within  the  scan  range  of  the  tunable  laser.  Fig.  3. 
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4.0  ComDarison  of  Code  Predictions  With 
Experiments 


Figure  3.  Schematic  of  the  arcjet  and  the  two  diagnostic 
laser  beams.  The  dot  at  the  intersection  of  the  laser 
beams  represents  the  detection  volume,  as  defined  by  the 
detection  optics. 

Since  the  distribution  function  in  arcjet  flows 
is  essentially  Maxwellian,  the  Doppler  line  profile  is 
Gaussian.  The  radial  component  of  the  mean  veloc¬ 
ity,  <Vy  >,  is  determined  directly  from  the  perpen¬ 
dicular  beam.  For  Gaussian  line  profiles  the  velocity 
is  determined  from  the  shift  of  the  measured  line  cen¬ 
ter  relative  to  the  reference  line  from  the  optogalvanic 
cell  using  Eq.  3.2.  The  axial  component,  <  Vx  >, 
is  then  determined  from  the  mean  velocity  measured 
along  direction  <j), 

<Vx>={<V6>  -  <Vy>  sm(j))/cos(j)  . 

(3.5) 

In  this  way  both  radial  and  axial  components  of  the 
mean  velocity  can  be  measured  simultaneously. 

Since  the  LIF  technique  measures  the  velocity 
distribution  function,  it  is  also  possible,  in  principle, 
to  determine  the  temperature  of  the  flow  from  the 
second  moments  of  the  measured  Doppler  line  pro¬ 
file.  In  practice,  we  found  considerable  differences 
in  the  temperatures  determined  in  this  way  from  the 
Doppler  profiles  obtained  at  different  angles,  result¬ 
ing  in  a  higher  temperature  in  the  direction  parallel  to 
the  arcjet  axis.  Since  the  second  moment  is  a  measure 
of  the  spread  in  velocity,  it  will  include  fluctuations  of 
the  mean  velocity,  in  addition  to  the  spread  in  ther¬ 
mal  velocity,  when  averaged  over  time,  as  i^ui^  by 
the  phase-lock  instrumentation.  Fluctuations  in  the 
mean  velocity  can  occur  as  a  result  of  unsteady  op¬ 
eration  of  the  arcjet  or  from  turbulence  in  the  flow. 
Thus,  temperature  measurements  obtained  from  the 
Doppler  profiles  represent  an  upper  limit  for  the  true 
temperature  in  the  flow. 


Predictions  of  the  UTSI  arcjet  simulation  code 
have  been  compared  with  experimental  results  for 
several  different  arcjets.  A  brief  overview  of  these  re¬ 
sults  will  be  given  here;  a  more  complete  description 
can  be  found  in  the  cited  references. 

4.1  Comparisons  of  Code  Predictions  with 
Experiments 

The  UTSI  equilibrium  arcjet  model,  which  as¬ 
sumes  local  thermodynamic  equilibrium  (LTE),  was 
first  used  to  simulate  a  30  kW  ammonia  arcjet  that 
was  tested  by  the  Jet  Propulsion  Laboratory  (JPL) 
[13].  The  numerical  model  overestimated  specific  im¬ 
pulse  (Isp)  by  about  12  percent,  and  the  numerical 
anode  to  cathode  voltage  drop  was  about  20  percent 
high.  A  major  conclusion  of  this  study  was  that  im¬ 
proved  procediures  for  calculating  transport  proper¬ 
ties  were  necessary  to  improve  the  accuracy  of  the 
model. 

Next,  a  second  version  of  the  UTSI  equilibrium 
arcjet  model  was  developed  for  hydrogen  propellant. 
The  hydrogen  and  ammonia  versions  of  the  arcjet 
model  were  then  used  to  compare  the  performance  of 
a  hydrogen  and  an  ammonia  arcjet  of  the  same  con¬ 
figuration  and  at  the  same  electrical  power  [21].  The 
30  kW  arcjet  configuration  tested  at  JPL  and  mod¬ 
eled  previously  was  chosen  for  this  work.  Several  im¬ 
provements  had  been  made  to  the  code  since  the  first 
simulation.  The  predicted  Isp  was  5.6  percent  lower 
than  the  earlier  calculation,  but  was  still  11  percent 
higher  than  the  arcjet  data.  As  expected,  when  the 
arcjet  configuration  was  modeled  with  hydrogen  pro¬ 
pellant,  the  predicted  Isp  was  50  percent  higher  than 
the  Isp  for  ammonia  (21).  The  effects  of  transport 
properties  on  arcjet  performance  were  also  investi¬ 
gated  in  this  work.  Numerical  studies  were  made  to 
assess  the  effect  of  transport  properties  on  the  pre¬ 
dicted  performance  of  the  hydrogen  arcjet.  These 
studies  indicated  that  diffusive  transport  is  very  sig¬ 
nificant  even  in  the  supersonic  part  of  the  flow,  and 
that  relatively  small  changes  in  transport  propertira 
can  have  a  significant  effect  on  performance.  This 
work  also  showed  that  non-equilibrium  recombination 
chemistry  can  have  a  large  effect  on  transport  proper¬ 
ties.  It  was  concluded  that  finite  rate  chemistry  cal¬ 
culations  are  necessary  if  accurate  arcget  performance 
is  to  be  calculated. 

Next,  the  UTSI  equilibrium  arcjet  model  was 
tested  for  three  different  nozzle  configurations  of  a 
high  power  hydrogen  arcjet  at  a  number  of  power 
levels  and  mass  flow  rates  [15].  The  three  nozzle  con¬ 
figurations  differed  in  constrictor  length,  constrictor 
diameter,  nozzle  half-angle  and  nozzle  expansion  ra¬ 
tio.  Experimental  data  for  these  configurations  were 
provided  by  NASA  Lewis  Research  Center.  The  gen¬ 
eral  trends  of  the  data  were  predicted  by  the  model. 
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However,  both  Isp  and  efficiency  were  overpredicted 
by  the  model,  the  Isp  by  as  much  as  10  percent  and 
the  efficiency  by  as  much  as  25  percent.  The  good 
qualitative  agreement  between  the  model  calculations 
and  the  data  indicated  that  the  general  procedures 
used  in  the  model  properly  accounted  for  the  essen¬ 
tial  physics  of  the  flow  and  heat  addition.  The  lack 
of  quantitative  agreement  indicated  that  the  numer¬ 
ical  values  used  for  transport  properties  may  have 
been  in  error,  and  that  one  or  more  loss  mechanisms 
were  being  neglected  or  underestimated.  These  re¬ 
sults  supported  the  need  to  include  both  a  finite  rate 
chemistry  model  to  account  for  frozen  flow  losses,  and 
a  two-temperature  model  for  electrons  and  the  heavy 


gas. 

After  multiplexed  LIF  data  on  a  small  water- 
cooled  argon  arcjet  was  collected  (11],  detailed  exit 
plane  profiles  of  axial  velocity,  radial  velocity,  and 
temperature  were  available  for  comparison  with  the 
arcjet  model.  An  argon  version  of  the  UTSI  equi¬ 
librium  arcjet  model  was  developed  for  this  work 
[25].  The  model  predicted  axial  velociti«  consis¬ 
tently  higher  than  those  determined  experimentally, 
although  the  calculated  and  experimental  velocity 
profiles  did  have  a  similar  shape.  The  calculated 
radial  velocity  matched  the  experiment  results  very 
well.  The  model  exit  temperatures  were  twice  those 
calculated  from  the  experiment  data.  The  differences 
between  the  experiment  and  the  model  results  were 
believed  to  be  the  result  of  underprediction  of  the  en¬ 
ergy  losses  by  the  UTSI  arcjet  model.  Energy  losses 
from  radiation  and  wall  conduction  in  a  small  water 
cooled  arcjet  should  be  a  significant  fraction  of  the 
input  power.  The  numerical  solution  showed  essen¬ 
tially  no  energy  loss.  Radiative  transport  was  not 
included  in  the  model  because  the  operating  pressure 
was  expected  to  be  too  low  for  the  radiation  of  argon 
to  be  significant.  It  was  concluded  that  radiation  en¬ 
ergy  transfer,  a  two-temperature  model  for  electrons 
and  the  heavy  gas,  and  finite  rate  chemistry  would  be 
necessary  to  better  predict  experimental  data. 

Next,  the  hydrogen  equilibrium  code  was  used 
to  simulate  some  interesting  arcjet  experiments  per¬ 
formed  at  The  Universitat  Stuttgart  Institut  fur 
Raumfahrt  Systeme  (IRS)  [8].  The  experiments 
utilized  two  arcjets;  one  with  a  high  temperature 
radiation-cooled  anode,  and  a  similar  arcjet  which 
used  a  water  cooled  anode.  Some  difficulty  was  expe¬ 
rienced  in  convergence  of  the  water-cooled  solution. 
Many  attempts  were  made  to  improve  conve^ence 
with  relatively  little  success.  The  best  solution  for  the 
water-cooled  arcjet  was  used  for  companrons  with  the 
solution  for  the  radiation-cooled  arcjet.  The  ex¬ 
periments  showed  that  the  efficiency  of  the  radiation- 
cooled  arcjet  was  41.4  percent  but  the  efficiency  of 
the  water-cooled  arcjet  was  only  18.5  percent  [8^  Our 
simulations  of  this  experiment  clearly  revealed  the  in¬ 
adequacies  of  the  equilibrium  model  and  the  ne>insric 
modification  of  the  electrical  conductivity.  The  effi¬ 
ciency  prediction  for  the  radiation-cooled  arqet  was 


45.1  percent,  in  reasonable  agreement  with  the  exper¬ 
iment,  but  the  prediction  for  the  water-cooled  arcjet 
was  32  percent,  nearly  a  factor  two  higher  than  the 
experiment  [8]. 

Recently,  the  UTSI  equilibrium  code  was  used 
to  simulate  a  1  kW,  ammonia  arcjet  whose  perfor¬ 
mance  was  measured  at  NASA  Lewis  Research  Cen¬ 
ter  [26|.  Thrust  experiments  were  performed  to  deter¬ 
mine  specific  impulse  for  a  range  of  electrical  power 
inputs  and  propellant  flows.  Computer  simulations 
predicted  specific  impulse  for  the  same  arcjet  geome¬ 
try  using  identical  power  inputs  and  propellant  fio^s. 
Figure  4  compares  the  specific  impulse  of  the  model 
with  the  experimental  measurements  as  a  function  of 
specific  energy.  Specific  energy  is  the  electrical  energy 
deposited  in  the  flow  per  unit  mass  of  propellant. 


Figure  4.  Comparison  of  measured  and  predicted  specific 
impulse  as  a  function  of  specific  energy  for  a  1  kW  am- 
monia  arcjet. 


Both  the  calculated  and  the  measured  specific 
impulse  appear  to  increase  linearly  with  specific  en¬ 
ergy.  However,  the  predicted  values  significantly  ^- 
ceed  the  experimental  values.  This  lea^  to  the 
conclusion  that  significant  energy  losses  have  been 
neglected.  The  assumption  of  thermal  and  chemi¬ 
cal  equilibrium  neglects  important  loss  mechanisms. 
Non-equilibrium  effects  in  the  nozzle  prevent  thermal 
energy  that  has  been  expended  in  dissociation  and 
ionization  from  being  converted  into  kinetic  energy 
during  expansion.  The  true  flow  becomes  frozen  in  a 
particular  chemical  state,  but  the  equilibrium  model 
allows  the  energy  to  be  recovered.  The  equilibrium 
model  also  does  not  account  for  heat  conduction  to 
the  wall  by  diffusion  of  high  energy  electrons  from 
the  high  temperature  core  through  the  cold  bound¬ 
ary  layer.  The  assumption  of  equilibnum  ako  alters 
the  character  of  the  flow  by  changing  the  el^tncal 
conductivity,  which  controk  the  distribution  of  power 
absorption. 
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Each  time  a  solution  from  the  UTSI  equilib¬ 
rium  code  was  compared  to  experimental  data,  it  was 
concluded  that  the  equilibrium  assumption  should  be 
relaxed.  A  new  non-equilibrium  version  of  the  hy¬ 
drogen  arcjet  code  was  developed  which  included  fi¬ 
nite  rate  chemistry  and  separate  energy  equations  for 
electrons  and  heavy  particles.  Predictions  from  the 
new  non-equilibrium  code  were  compared  with  a  high 
powered  arcjet  experiment  run  at  NASA  Lewis  Re¬ 
search  Center  [12].  Some  difficulty  was  experienced 
with  the  convergence  of  this  solution,  but  the  predic¬ 
tions  of  specific  impulse  and  efficiency  were  in  rea¬ 
sonable  agreement  with  the  experimental  thrust  and 
efficiency.  The  new  non-equilibrium  results  under¬ 
predicted  the  experimental  results  by  as  much  as  the 
equilibrium  code  overpredicted  them. 

5.0  Multiplexed  LIF  Experiments 

The  multiplexed  LIF  technique  was  first  used  at 
UTSI  to  probe  the  plume  of  a  300  Watt  argon  arcjet 
[11].  The  purpose  was  to  demonstrate  the  method 
of  using  two  multiplexed  laser  beams  to  measure  two 
velocity  components  simultaneously  and  to  demon¬ 
strate  the  method  of  using  the  optogalvanic  effect  in 
a  hollow  cathode  discharge  as  a  stationary  reference. 
Argon  propellant  was  used  while  proving  the  con¬ 
cept  because  of  the  large  population  of  metastable 
excited  states  available,  but  the  LIF  technique  was 
later  applied  to  more  applicable  propellant  species. 
LIF  measurements  of  axial  velocity,  radial  velocity, 
and  temperature  were  obtained.  A  radial  scan  was 
made  across  the  nozzle  exit  one  millimeter  down¬ 
stream  from  the  exit  plane.  An  axiail  scan  was  ob¬ 
tained  downstream  of  the  nozzle  along  the  arcjet  cen¬ 
terline.  The  data  indicated  a  radial  boundary  layer 
and  a  normal  shock  downstream  from  the  nozzle.  The 
location  of  these  effects  was  consistent  with  the  barrel 
shock  and  Mach  disk  visually  observed  in  the  plume. 

Once  two-beam  multiplexed  LIF  measurements 
had  been  fully  demonstrated,  the  technique  was  mod¬ 
ified  to  probe  the  plume  of  a  one  kilowatt  arcjet  us¬ 
ing  simulated  ammonia  propellant.  In  preparation 
for  the  LIF  measurements,  several  experiments  were 
performed  on  the  plume  of  the  ammonia  arcjet  using 
emission  spectroscopy  [29].  The  temperature  of  the 
plume  at  the  exit  plane  was  determined  with  three  in¬ 
dependent  methods:  1)  Doppler  width  of  the  hydro¬ 
gen  Balmer  alpha  line,  2)  intensity  ratios  of  the  hy¬ 
drogen  B aimer  alpha  and  beta  lines,  and  3)  rotational 
line  intensities  of  the  NH  radical.  The  three  temper¬ 
ature  measurements  were  in  fairly  close  agreement, 
but  the  spatial  trends  did  not  agree.  Emission  spec¬ 
troscopy  was  also  used  to  make  measurements  of  the 
axial  velocity  of  the  plume  by  observing  the  Doppler 
shift  of  the  hydrogen  Balmer  alpha  line.  Rapid  accel¬ 
eration  of  the  plume  was  observed  downstream  of  the 
exit  plane,  indicating  an  underexpanded  flow.  The 
emission  velocity  data  is  somewhat  difficult  to  inter¬ 
pret  because  emission  spectroscopy  is  a  line-of-sight 


measurement,  but  it  is  believed  that  emission  exper¬ 
iments  can  be  used  to  make  an  estimation  of  axial 
plume  velocity  when  LIF  is  either  impractical  or  not 
available. 

The  emission  results  were  later  compared  to  LIF 
measurements  performed  at  identical  conditions  [30] . 
The  LIF  experiments  were  performed  by  exciting  the 
hydrogen  Balmer  alpha  transition.  The  LIF  measure¬ 
ments  confirmed  the  acceleration  of  propellant  down¬ 
stream  of  the  nozzle  exit  as  observed  in  emission.  Sev¬ 
eral  experimental  factors  reduced  the  quality  of  the 
hydrogen  LIF  data  below  that  obtained  with  argon. 
The  width  of  the  hydrogen  Doppler  profile  is  similar 
to  the  maximum  scan  width  of  the  laser  (30  GHz) 
which  makes  the  profile  difficult  to  curve  fit.  In  ad¬ 
dition,  the  signal-to-noise  ratio  was  high  because:  1) 
the  LIF  signal  obtained  from  the  Balmer  alpha  tran¬ 
sition  was  much  less  than  the  signal  acquired  from  the 
argon  transition,  and  2)  there  was  a  large  amount  of 
background  laser  scatter. 

Additional  LIF  measurements  in  the  ammonia 
ELTcjet  plume  were  made  to  obtain  data  for  compari¬ 
son  with  the  UTSI  arcjet  simulation  code  [26].  The 
signal-to-noise  ratio  was  greatly  improved  by  using 
more  sensitive  optical  detection  and  by  reducing  laser 
scatter.  The  LIF  technique  was  used  to  measure 
the  plume  velocities  and  temperatures  in  radial  scans 
across  the  exit  plane.  In  these  experiments,  both 
the  hydrogen  Balmer  alpha  line  and  a  transition  of 
atomic  nitrogen  were  used  in  the  multiplexed  LIF 
measurements.  The  nitrogen  profiles  are  much  nar¬ 
rower  than  the  scan  width  of  the  laser,  and  there  were 
no  problems  fitting  the  data.  In  addition,  the  fluores¬ 
cence  could  be  detected  at  a  different  wavelength  than 
the  laser,  virtually  eliminating  the  effects  of  scattered 
laser  radiation.  The  comparisons  between  the  equi¬ 
librium  model  and  the  data  are  shown  in  Figs.  5-7. 

An  interesting  result  of  these  experiments  was 
that  atomic  hydrogen  and  atomic  nitrogen  in  the  flow 
had  different  velocities  (Figs.  5  and  6). 

The  hydrogen  is  expanding  radially  much  faster 
than  the  nitrogen,  but  the  nitrogen  appears  to  have 
a  slightly  larger  axial  component  of  velocity.  It  is 
believed  that  the  hydrogen  is  diffusing  outward  more 
rapidly  because  it  is  lighter  than  nitrogen. 

In  general,  the  temperatures  measured  with  LIF 
are  less  accurate  than  the  velocity  measurements.  It 
is  more  difficult  for  the  fitting  program  to  determine 
accurate  widths  than  accurate  profile  centers  from  the 
experimental  data.  The  width  of  the  profile  obtained 
from  the  fitting  program  is  sensitive  to  changes  in 
amplitude  and  baseline  offset  of  the  LIF  signal,  while 
the  determination  of  the  center  of  the  profile  is  rela¬ 
tively  insensitive  to  these  effects.  Useful  temperature 
data  was  not  obtained  from  atomic  hydrogen  because 
of  the  extreme  width  of  the  profile.  A  radial  scan  of 
temperatures  determined  using  the  atomic  nitrogen 
data  is  shown  in  Fig.  7. 
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Figure  5.  Measured  axial  velocities  in  the  ammonia  arcjet 
plume  one  millimeter  from  the  nozzle  exit.  LIF  measure¬ 
ments  are  shown  for  both  atomic  hydrogen  and  atomic 
nitrogen. 


Figure  6.  Measured  radial  velocities  in  the  ammonia  ar¬ 
cjet  plume  one  millimeter  from  the  nozzle  exit.  LIF 
measurements  are  shown  for  both  atomic  hydrogen  and 
atomic  nitrogen. 


Figure  7.  Measured  temperatures  obtained  from  the  ni¬ 
trogen  LIF  experiments  one  millimeter  from  the  nozzle 
exit.  Data  are  shown  for  both  the  perpendicular  and  the 
angled  beams. 


The  temperatures  acquired  from  the  angled 
beam  are  consistently  higher  than  the  temperatures 
acquired  from  the  perpendicular  beam.  This  effect  is 
believed  to  be  caused  by  broadening  of  the  Doppler 
profile  due  to  spatial  velocity  gradients  in  the  flu¬ 
orescence  collection  volume.  Previous  experiments 
performed  on  this  nozzle  indicate  that  the  propellant 
undergoes  rapid  axial  acceleration  after  it  leaves  the 
nozzle.  The  velocity  gradients  should  be  larger  in  the 
axial  direction  than  in  the  radial  direction.  Therefore, 
the  angled  beam  will  sample  a  larger  variation  in  ve¬ 
locity  than  the  perpendicular  beam,  and  the  Doppler 
profile  detected  by  the  angled  beam  will  be  broadened 
more  than  the  profile  of  the  perpendicular  beam. 

6.0  Recent  Calculations  Using 
the  Non-eouilibrium  Code 

Several  significant  changes  have  been  made  to 
the  non-equilibrium  code.  In  the  present  version  of 
the  code,  only  the  value  of  total  current  is  fixed,  and 
the  total  power  input  is  calculated  by  the  code.  This 
is  now  possible  because  the  distribution  of  electrical 
conductivity  within  the  flow  calculated  for  the  non¬ 
equilibrium  plasma  is  more  accurate  than  that  for 
the  LTE  case.  The  assumption  that  the  current  flow¬ 
ing  into  the  anode  decreased  linearly  with  axial  posi¬ 
tion  has  also  been  dropped.  Now,  the  current  is  only 
forced  to  flow  into  the  anode  in  a  normal  direction, 
consistent  with  a  constant  anode  potential,  and  the 
distribution  of  current  into  the  anode  is  calculated 
by  the  code.  These  changes  result  in  a  more  realis¬ 
tic  simulation  and  have  also  improved  convergence, 
although  the  number  of  iterations  required  to  reach 
convergence  has  increased. 

Calculations  were  made  using  the  new  version  of 
the  non-equilibrium  code  to  compare  with  two  arcjet 
experiments  reported  by  docker  and  Auweter-Kurtz 
at  the  Institut  fur  Raumfahrt  Systeme  in  Stuttgart. 
These  are  the  radiation-cooled  MARC  thruster  [31] 
and  the  water-cooled  TTl  thruster  [32].  Moeller  has 
compared  the  UTSI  equilibrium  arcjet  model  with 
these  experiments  [8].  These  experiments  were  se¬ 
lected  for  comparison  because  they  have  similar  con¬ 
figurations  and  operating  parameters,  but  exhibit  a 
large  difference  in  performance.  Also,  Miller  and 
Martinez-Sanchez  have  used  these  experiments  for 
comparison  with  their  arcjet  calculations  [33].  The 
geometries  of  these  thrusters  are  given  in  Table  1, 
and  the  computational  grid  used  for  both  thrusters  is 
shown  in  Fig.  8.  Convergence  is  a  necessary,  but  not 
sufficient,  requirement  for  a  satisfactory  simulation. 
Various  tests  have  been  made  to  evaluate  the  qual¬ 
ity  of  the  solution.  Several  well-converged  solutions 
have  been  obtained  with  the  non-equilibrium  arcjet 
model.  Some  examples  of  these  tests  are  shown  for  a 
well-converged  solution  for  the  MARC  thruster  using 
the  “fast”  set  of  kinetic  rates,  to  be  described  later. 
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Table  1.  Geometric  parameters 


for  the  MARC  and  TTl  hydrogen  arcjet  thrusters. 


Constrictor  Length  (mm) 
Constrictor  Diameter  (mm) 
C^hode  Tip  Half  Angle  (deg) 
Nozzle  Area  Ratio 
Xhode/Cathode  Gap  (mm) 
Nozzle  Half  Angle  (deg) 


MARC  Thruster 

5.0 

2.5 
30 
100. 

2. 

17.5 


TTl  Thruster 

5.0 

2.5 

30 

100. 

2. 

Contoured 


Figure  8.  Computational  grid  used  for  both  the  MARC 
and  TTl  hydrogen  arcjet  thrusters. 

One  test  of  the  solution  is  to  determine  how  well 
mass  is  being  conserved.  The  radial  integral  of  mass 
flow  from  the  converged  solution  is  shown  in  Fig.  9  as 
a  function  of  axial  location.  The  mass  flow  is  typically 
constant  to  within  a  few  percent. 


Figure  9.  Radial  integral  of  the  calculated  mass  flow  as  a 
function  of  axial  position  for  a  well-converged  solution. 
The  correct  mass  flow  is  0.1  g/s. 


Power  balances  are  checked  by  comparing  the 
total  energy  flux  in  the  flow  (enthalpy  plus  kinetic 
energy)  with  the  cumulative  electrical  powCT  input 
and  the  cumulative  wall  losses.  These  quantities  are 
shown  in  Fig.  10  together  with  the  integrated  error. 
The  error  shown  in  this  figure  is  the  sum  of  the  resid¬ 
ual  errors  at  each  grid  point  upstream  of  the  axial 
position.  The  residual  errors  are  clearly  small  com¬ 
pared  with  the  energy  input. 

The  quality  of  the  solution  for  the  electromag¬ 
netic  field  is  evaluated  by  examining  voltage  and 
current  distributions.  The  current  distributions  are 
shown  in  Fig.  11. 


Figure  10.  Comparison  of  the  energy  flux  in  the  flow  with 
electrical  power  input,  wall  power  loss,  and  the  residual 
error. 


Figure  11.  Current  distributions.  The  radial  integral  of 
the  axial  component  of  current  and  the  cumulative  ar¬ 
rant  into  the  anode  sum  to  the  sp^ified  80  A  total.  Also 
shown  is  the  current  distribution  in  the  anode. 
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The  sum  of  the  radial  Integral  of  the  axial  com¬ 
ponent  of  the  current  and  the  cumulative  current 
flowing  into  the  anode  is  essentially  constant  at  the 
prescribed  value  of  80  amps.  The  current  distribution 
into  the  anode  is  strongly  peaked  near  the  cojistnctor 
exit,  but  some  anode  current  extends  over  the  length 

of  the  nozzle.  ry  v  p 

For  the  steady-state  case  considered,  V  X  li 
must  equal  zero  (Eq.  2.7).  This  means  that  the  line 
integral  of  electric  field  around  any  closed  path  must 
vanfeh.  However,  the  calculated  cathode  to  anode 
voltage  is  not  completely  independent  of  the  path 
along  which  the  voltage  is  calculated.  How  well  this 
requirement  is  preserved  by  the  code  can  be  evalu¬ 
ated  by  comparing  the  voltage  between  the  anode  and 
the  centerline  as  calculated  along  two  different  paths. 
One  path  begins  at  a  point  on  the  anode,  then  goM 
one  cell  toward  the  inlet  and  then  radially  inward  to 
the  axis.  The  other  path  goes  from  the  point  on  the 
anode  radially  to  the  axis  and  then  one  cell  toward 
the  inlet.  The  results  are  shown  in  Fig.  12. 

There  is  good  agreement  for  the  two  paths  ex¬ 
cept  near  the  constrictor/nozzle  comer  at  an  axial 
location  of  0.0155  m.  The  cathode-to-anode  voltage, 
like  the  total  electrical  power  input,  is  calculated  by 
the  code.  This  voltage  calculated  from  a  line  Integra 
from  the  cathode  tip  along  the  centerline  to  the  axial 
position  of  the  constrictor/nozzle  comer,  then  radi¬ 
ally  to  the  anode  and  then  along  the  anode  is  some¬ 
what  different  than  if  the  calculation  path  is  abng 
the  centerline  and  then  radially  to  the  anode.  This 

effect  is  shown  in  Fig.  13.  ic 

In  both  cases  the  calculated  anode  voltage  is 

essentially  constant  along  the  nozzle,  and  the  current 
enters  the  anode  normal  to  the  wall. 


■  -  PcHn  loword  Wet  on  anode,  then  lo  centerline 
#  -  Path:  to  cenlefSne  from  anode,  then  toword  Wet 


♦  -  A»«ol  on  centerfine.  then  rorfd  to  onode 
■  -  Axtd  on  center  lo  throat,  rodd  lo  onode,  then  a«ai 
^  -  Experimenl 


- 1^5  .02  .025  Ti  .035  .04  .045  .05 

Axid  ds  lance 

Figure  12.  Anode  to  axis  voltage  calculated  from  the  line 
integral  along  different  paths. 


.01  .015  .02  .025  .03  .035  .04*  .045  .05 

Axial  distance  -m 

Figure  13.  Voltage  from  the  cathode  tip  to  a  point  on 
the  anode  surface  as  calculated  along  two  different  paths. 
The  experimental  cathode-to-anode  voltage  of  112  V  is 
shown  for  comparison. 

fi1  MARC  Thruster  Simulations 

Two  calculations  were  made  for  the  MARC 
thruster  using  different  rates  for  the  ion^tion  of  H 
atoms  by  electrons  (Table  2).  The  “=1°^ 
in  the  model  was  presented  at  the  23rd  lEPC  [  ] 

with  the  recombination  rate  slightly  modified  to  fit 
the  data  from  Park  [22].  The  “fast”  rate  us^  t^ 
reverse  rate  of  the  recombination  reaction  obtain^ 
through  the  equlibrium  constant  evaluated  at  the 
electron  temperature.  The  recombination  reaction 
rate  was  taken  from  McCay  (341  as  were  the  rates 
for  the  other  reactions  in  both  sets. 

The  calculated  performance  for  the  fa^  and 
“slow”  rates  and  for  an  equlibrium  calculation  are 
compared  with  the  experimental  data  m  Table  i. 

The  calculation  with  the  fast  chemical  rates 
most  closely  matched  the  data.  The  thrust  is  about 
6  5  percent  low  and  the  calculated  power  and  volt¬ 
age  about  9  percent  high.  The  equilibrium  calcula¬ 
tion  overpredicts  thrust  by  more  than  25  percent  Md 
power  by  more  than  50  percent.  The  agr^rnent  be¬ 
tween  the  power  calculated  from  the  product  of  cur- 
rlTand  voltage  and  that  calculated  from  the  volume 
integral  of  is  also  better  for  the  case  run  with 

the  “fast”  chemistry.  . 

The  “slow”  chemistry  calculation  gives  about 

4  percent  less  thrust  and  dissipates  about  8  Percent 
more  power  than  the  “fast”  cheinistry  case.  This  rel¬ 
atively  small  difference  is  somewhat  suprwing  consid¬ 
ering  the  large  differences  in  predicted  electron  con¬ 
centration  shown  in  Fig.  14. 
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Table  2.  Fast  and  slow  reaction  rates  [12  and  34]. 


Forward  rate  =  AT  ^  exp(£^/T) 


Reaction 

A 

n 

E 

1  H  +  H-\-M  =  H2  +  M 

6.40e  +  17 

1.0 

0. 

2  H+  +  e  +  M  =  H  +  M 

5.26e  +  26 

2.5 

0. 

3  H2  +  6  — ^  H  H  -h  6 

1.91e  +  ll 

-1.0 

-203. 

4  H+  +  e-^H  +  hu 

3.77e  + 13 

0.58 

0. 

fast  rates 

5  e+  e  =  H e 

7.08e  +  39 

4.5 

0. 

slow  rates 

5r  +  e  +  e  ^  H  +  e 

2.19e  +  41 

5.06 

0. 

5i  ff  +  e  — »•  +  e  +  e 

1.81e+14 

0.0 

-244750, 

M  any  third  body 
— ►  one  way  reaction  -  no  reverse  rate 
=  reversible  reaction 

-  reverse  rate  obtained  from  equilibrium  constant 
5r  recombination 
5i  ionization 

Reaction  1  and  2,  t  =  gas  temperature 
Reaction  3,  4,  and  5,  T  =  electron  temperature 


Table  3.  Comparison  of  calculated  results  from  the  non-equilibrium  code  using  two  different  sets  of 
reaction  rates  with  an  equilibrium  calculation  and  experimental  data  for  the  MARC  thruster. 


Area  (m^) 
Mass  (kg) 
Pinlet  (Atm.) 
Thrust  (N) 
isp  (sec) 
ke  (W) 
Efficiency 
1(A) 

E(V) 

El  (kW) 

J^  /a  (kW) 


Fast  Rate 

.673  X  10"^ 
.100  X  10-3 
1.01 
0.828 
845.3 
3431. 
0.362 

80.2 
125. 
10.0 
9.863 


Calculation 
Slow  Rate 
.673  X  lO-'^ 
.100  X  10-3 
1.00 
0.795 
811.3 
3161. 
0.333 

80.2 
140. 

11.2 
10.68 


Equilibrium 

.673  X  10-^ 
.100  X  10-3 
1.08 
1.103 
1125.7 
6086. 
0.642 

80.2 
170. 
13.6 
14.10 


Data 

.673  X  10-^ 
.100  X  10-3 

0.866 

882.8 

0.414 

80.2 

112.9 

9.053 
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Figure  14.  Calculated  electron  concentrations.  ^  The 
“fest”  chemistry  results  are  shown  in  (a)  and  the  slow 
chemistry  results  are  shown  in  (b). 


However,  the  electron  temperature  increases 
when  the  electron  concentration  decreases,  as  shown 

in  of  the  inverse  relation  between  elec¬ 

tron  concentration  and  electron  temperature  is  a 
smaller  change  in  electrical  conductivity  than  would 
or^if  only  the  electron  concentration  or  tempera- 

ture  changed. 


Figure  15.  Calculated  electron  temperatures.  The  f^t 
chemistry  results  are  shown  in  (a)  and  the  slow  chem¬ 
istry  results  are  shown  in  (b). 


9.  TTl  Thruster  Simulations 

A  converged  solution  for  the  TTl  thruster  wm 
btained  for  the  “slow”  chemistry  which  agrees  quite 
rell  with  the  experimental  data  for  this  case  (Ta 

The  configuration  of  both  these  throsters  are 
he  same,  except  that  the  TTl  thruster  has  a  con- 
oured  nozzle.  The  conical  MARC  nozzle  was  used 


Table  4.  Comparison  of  calculated  results  from  the 

rates  with  experimental  data  for  the  TTl 


Area 

Mass  (kg) 
Pinlet  (Atm.) 
Thrust  (N) 
isp  (sec) 
ke  (kW) 
EfRciency 
1(A) 

E(V) 

El  (kW) 
JVo-  (kW) 


Calculation  Data 

.673  X 

.673  X  10"" 

.100  X  10~® 

.100  X  10"' 

0.68 

— 

0.603 

0.626 

615.6 

638.2 

3431. 

— 

0.173 

0.185 

61.3 

61.3 

160. 

160. 

9.8 

9.814 

9.404 

— 
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for  this  simulation  because  the  new  algorithm  used 
to  obtain  the  zero  normal  gradient  boundary  condi¬ 
tion  for  the  magnetic  field  on  the  anode  has  not  yet 
been  generalized  to  work  for  a  contoured  nozzle.  The 
“fast”  chemistry  and  the  equlibrium  cases  failed  to 
converge  for  the  TTl  case,  with  the  residual  errors 
reaching  a  stable  level  with  a  rather  regular  racilla- 
tion  about  this  level.  The  failure  of  the  equlibrium 
case  is  rather  surprising,  since  the  terms  which  seem 
most  likely  to  introduce  the  instabilities  do  not  occur 
in  this  case. 

fi-.l  Status  oF  the  Computational  Code 

The  UTSI  model  has  been  able  to  provide  con¬ 
verged  solutions  to  four  out  of  six  cases  attempted. 
Convergance  is  slow  and  the  residual  errors  oscillate. 

It  has  taken  up  to  60,000  iterations  to  obtain  a  satis¬ 
factory  solution.  This  required  about  50  hours  on  a 
6  megaflop  IBM  320  RISC  6000  machine. 

At  this  point  we  feel  that  the  algorithm  is  ^v- 
ing  solutions  which  are  trustworthy  enough  to  make 
some  evaluations  of  the  physics  included  in  the  model. 
First,  we  neglect  the  excitation  energy  of  the  hydro¬ 
gen  atom.  Park  [22]  proposes  the  assumption  that 
the  electronic  excitation  levels  are  in  equilibrium  with 
the  electron  temperature.  Adopting  this  assumption 
would  have  the  effect  of  increasing  the  effective  elec¬ 
tron  specific  heat  and  reducing  the  sensitivity  of  the 
electron  temperature  to  variations  in  energy  input. 
This  should  make  the  equations  more  stable  and  pro¬ 
mote  faster  convergence.  The  effect  of  this  change  on 
the  ability  of  the  model  to  predict  experimental  data 
is  unknown  at  this  time.  Second,  the  ability  to  pre¬ 
dict  the  correct  voltage  and  power  depends  on  having 
the  correct  spatial  distribution  of  the  electrical  con¬ 
ductivity.  There  appear  to  be  significant  regions  in 
the  flow  where  the  gas  deviates  from  an  ideal  plasma, 
as  determined  by  the  value  of  log  A  [17].  This  occurs 
where  the  electron  temperature  is  low  and  the  elec¬ 
tron  density  is  relatively  high.  The  points  in  the  con¬ 
verged  solution  for  the  MARC  thruster  wth  “fast’ 
chemistry  where  log  A  <  4  is  shown  in  Fig.  16. 

There  is  a  need  to  develop  a  proper  and  efficient 
method  to  calculate  electrical  and  thermal  conduc¬ 
tivities  in  a  non-ideal,  non-equilibrium  plasma.  The 
electron-ion  collision  cross-section,  which  is  a  major 
factor  in  the  electrical  conductivity,  also  appears  in 
other  transport  terms  and  in  the  rate  at  which  elec¬ 
trons  exchange  energy  with  the  heavy  spwies.  Modi¬ 
fications  to  account  for  plasma  non-ideality  will  have 
consequences  that  extend  beyond  the  effects  of  el^- 
trical  conductivity.  The  relative  insensitivity  of  the 
results  to  the  chemical  rates  needs  to  be  studied^  in 
greater  detail  to  determine  if  it  is  a  general  conclusion 
or  one  specific  to  the  cases  tested. 


Figure  16.  Points  in  the  converged  solution  for  the 
MARC  thruster  with  “fast”  chemistry  where  log  A  <  4. 

In  principal,  the  UTSI  non-equlibrium  arcjet 
model  could  be  applied  to  arcjets  using  the  hydrogen- 
nitrogen  propellants  (ammonia  or  hydrazine)  after  a 
transport  and  reaction  model  for  this  system  is  de¬ 
veloped.  In  practice,  with  the  more  complex  system, 
convergance  problems  could  become  worse  and  the  ^- 
sumption  of  a  single  diffusion  coefficient  for  the  mix¬ 
ture  becomes  questionable.  If  diffusive  separation  of 
H  and  N  occurs  within  the  nozzle,  it  can  only  be  pre¬ 
dicted  if  indivdual  species  diffusion  coefficients  are 
obtained  and  the  equations  which  use  these  coeffi¬ 
cients  are  reformulated. 
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